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First-principles calculations of the epitaxial-strain phase diagram of perovskite SrCoOa are pre- 
sented. Through combination of the large spin-phonon coupling with polarization-strain coupling 
and coupling of the band gap to the polar distortion, both tensile and compressive epitaxial strain 
are seen to drive the bulk ferromagnetic-metallic (FM-M) phase to antiferromagnetic-insulating- 
ferroelectric (AFM-I-FE) phases, the latter having unusually low elastic energy. At these coupled 
magnetic-ferroelectric metal-insulator phase boundaries, cross responses to applied electric and mag- 
netic fields and stresses are expected. In particular, a magnetic field or compressive uniaxial stress 
applied to the AFM phases could induce an insulator-metal transition, and an electric field applied 
to the FM-M phase could induce a metal-insulator transition. 

PACS numbers: 71.30.-|-h, 77.80.-e, 75.85.-|-t, 63.20.-e 



Recent dramatic advances in the synthesis of complex 
oxides have opened exciting avenues for the theoreti- 
cal design and laboratory realization of new classes of 
functional oxide materials With tuning of a control 
parameter such as composition, artificial structuring or 
epitaxial strain, the competition of structural, magnetic 
and electronic degrees of freedom characteristic of com- 
plex oxides can result in two phases very close in energy 
but with very different structure and magnetic ordering. 
Desirable functional behavior can be obtained with the 
application of applied fields and/or stresses that produce 
switching between the two phases, a metal- insulator tran- 
sition triggered by electric or magnetic fields or ap- 
plied stress being of particular current interest. 0, 0] 

As an extension of previous work on epitaxial-strain- 
induced multifcrroicity in EuTiOa 0] and SrMnOs 
0], one strategy to identify the parent phase for a cou- 
pled magnetic-ferroelectric metal-insulator transition is 
to look for a high-symmetry metallic phase with a lat- 
tice distortion that is strongly destabilized by a change 
in magnetic order and, moreover, opens a band gap. This 
would result in a low-energy alternative phase, with dis- 
tinct magnetic order, a lattice distortion and a band 
gap, that could be stabilized by tuning of an appro- 
priate control parameter, establishing a phase boundary 
at which switching between the metallic and insulating 
phases would be possible. 

Here, we draw on the results of a first-principles sur- 
vey of the spin-phonon coupling in selected sets of mag- 
netic perovskite oxides [1] to identify SrCoOa as a suit- 
able parent phase. These calculations showed that cu- 
bic SrCoOs has a large spin-phonon coupling, with the 
lowest-frequency polar phonon being strongly destabi- 
lized by a change in magnetic order from ferromagnetic 
(FM) to the higher-energy C?-AFM ordering. In experi- 
mental observations, bulk SrCoOa is found to be cubic, 
with ferromagnetic order below Tc=280--305 K ^E^. 
Previously reported experiments on SrCoOa show that 
oxygen vacancies or chemical doping can lead to dra- 



matic chariges in structure, magnetic ordering and resis- 
tivity 12-13 • The large spin-phonon coupling seen 
in first-principles results suggests that similarly dramatic 
coupled structural, magnetic and metal-insulator phase 
transitions, and associated functional properties, could 
occur with epitaxial strain or isoelectronic substitution. 

In this paper, we investigate the ground-state epitax- 
ial strain phase diagram of SrCoOa from first principles 
and find phase transitions from the bulk FM-M phase to 
insulating AFM-FE phases occurring at accessible values 
of compressive and tensile epitaxial strain; furthermore, 
the energies of the strained phases are exceptionally low 
with respect to the bulk, promoting their experimental 
realization. Wc determine the magnetic ordering, the 
metallic (M) or insulating (I) character, and the polar- 
ization of ferroelectric (FE) phases as a function of epi- 
taxial strain. We find that it is the polar distortion that 
is responsible for opening the band gap in the insulat- 
ing phases, which offers the possibility of epitaxial-strain 
and electric-field tuning of the band gap in the AFM-FE 
phase. We predict large mixed magnetic-electric-elastic 
responses in the vicinity of the phase boundaries resulting 
from switching between the two phases, including metal- 
insulator transitions triggered by electric and magnetic 
fields and uniaxial stress. 

First-principles calculations were performed using 
density-functional theory within tl ie g eneralized gradient 
approximation GGA+U method J20| with the Perdew- 
Becke-Erzenhof parametrization [2l| as implemented in 
the Vienna Ab Initio Simulation Package (VASP- 
4.6) [22, 23]. The projector augmented wave (PAW) po- 
tentials [2^1 explicitly include 10 valence electrons for Sr 
{As^Ap^5s^), 9 for Co {M^4s^), and 6 for oxygen (25^2^4). 
The wave functions are expanded in a plane waves basis 
with an energy cut-off of 500 eV. We use the Dudarev (25| 
implementation with on-site Coulomb interaction U=2.5 
eV and on-site exchange interaction J//=1.0 eV to treat 
the localized d electron states in Co. The value of U 
chosen gives a value of 2.6/iB/f.u. and a lattice constant 
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TABLE I: Selected calculated phonon frequencies (cm^^) and 
relative total energy (meV/f.u.) of cubic SrCoOs at the cal- 
culated equilibrium lattice constants with FM (ao=3.842 A) 
and G-AFM (ao=3.885 A) orderings. The classification of the 
three zone-center mode eigenvectors as Slater, Last and Axe 
follows standard terminology for perovskites [321 ]. 
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ao= 3.842 Afor the cubic FM structure in good agree- 
ment with a previous calculation [26L 27| and recent mea- 

=3.835 A[9|) on stoichiometric 



surements (2.5^_b [U 



SrCoOa samples. This magnetic moment corresponds to 
an intermediate (^29623) spin state. 

The phonon frequencies of cubic SrCoOa with G- 
AFM and FM orderings were computed using the frozen 
phonon method at the F, M and R points in the irre- 
ducible Brillouin zone of the primitive pcrovskite unit 
cell, with atomic displacements of O.OOlaQ. To find the 
minimum-energy configuration in a given space group de- 
termined by freezing in one or more unstable modes of the 
cubic reference structure, we moved the atoms accord- 
ing to the conjugate-gradient algorithm until the resid- 
ual HcUman-Feynman forces were less than l.OmeV/A. 
An 8x8x8 Monkhorst-Pack k-point grid was used for 
5-atom unit cells, while a 4 x 4 x 4 grid was used for 
the \/2 X \/2 X \/2 (10-atom) supercell and \/2 x -\/2 x 2 
(20-atom) supercell To study the effects of epitaxial 



strain, we performed "strained-bulk" calculations l2£ 
in which calculations were performed for the periodic 
crystal with appropriate epitaxial constraints imposed on 
the in-plane lattice parameters, with all atomic positions 
and the out-of-plane lattice constant optimized. Epitax- 
ial strain is here defined relative to the computed lattice 
constant for the FM cubic perovskite structure (3.842A). 
Ferroelectric polarizations for the relaxed structures at 
each strain were computed by the Berry-phase method 



Selected phonon frequencies of cubic FM and G-AFM 
SrCoOa are given in Table HI A dramatic change of 
phonon frequency with magnetic ordering is seen for the 
Slater mode, which is stable in the ground state FM cubic 
structure but strongly unstable with the higher-energy 
G-AFM ordering. The corresponding eigenvectors of the 
Slater mode are (Sr,Co,0|| ,0_l) = (0.02,0.38,-0.32,-0.61) 
for FM and (0.00,0.43,-0.39,-0.58) for G-AFM, showing 
displacement of the i3-site Co cation relative to a fairly 
rigid oxygen octahedron network. The Axe mode is much 
less sensitive to the magnetic ordering, while the change 
in the Last mode is negligible. In addition, the oxygen 
octahedron rotation modes at R and M soften consider- 
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FIG. 1: (Color online) GGA-ff/ total energies of various 
structures labeled by space groups generated by freezing in 
selected modes of the primitive perovskite structure j33 |. 
as listed below. Calculations were performed at every 1% 
strain increment and interpolated. The energies of struc- 
tures with FM and G-AFM ordering are shown with red 
and blue symbols, respectively; the single FM curve is also 
labeled. Vertical black arrows at -0.72% and +2.04% strain 
indicate phase boundaries separating the Pna2i AFM-I-FE2 
phase, the PA/mmm FM-M phase, and the Pna2i AFM- 
I-FE^ phase, where FEz and FE^y denote the direction of 
ferroelectric polarization. PAmm — F J [001] , Amm2 = 



FjJllO], PA/mbm = Mj+JOOl], P46m = M}[QQ1] + FjJOOl], 

7i+[110] -HM3+JOOI], 



P2i/m = JOU] + A/+J100], Pbnm - ..4 
Pmc2i = i?J[110] +A/3+J001] +rj[110], Pna2i 
Mg+JOOl] -hFjJOOl]. 



Pj[110] + 



ably in the G-AFM state. 

For FM ordering, we find that the cubic structure is 
stable not only against the distortions described by single 
modes in Table HJ but also against the distortion into 
the orthorhombic Pbnm structure, common in perovskite 
oxides, generated by the combination i?^[110]-l-Af;^[001]. 
This is consistent with the experimental observation that 



bulk SrCoOa is cubic. 

In contrast, for G-AFM ordering, the cubic structure 
is unstable against distortion. In fact, the energy 
gains associated with several different distortions are 
comparable to 230meV/fu, the energy of the cu- 
bic G-AFM phase (ao=3.885 A) relative to the FM 
ground state phase. (P6nm=-|-89meV/fu, P4mm=-f 34, 
Amm2=+15, Pna2i{R+ [110] + M3+r4 [001])=+18, 
Pmc2i(i?|[110] M3+[001] + [110])=+14) 

In Fig. [U the epitaxial strain dependences of the total 
energies of the nonpolar PA/mmm structure and struc- 
tures with various polar, Jahn- Teller (JT), and rota- 
tional distortions, are shown for FM and G-AFM or- 
dering for strains from -1.9% (compressive) to -1-3.3% 
(tensile) strain. At 0% strain, the ground state cubic 
FM phase is the lowest energy state, while the G-AFM 
PA/mmm phase is much higher in energy. However, the 
polar instability of the latter phase (see Table |l| leads 
to large energy gains for polar G-AFM phases. At 0% 
strain, there is a gain relative to nonpolar PA/mmm of 
190 meV/f.u. with r4 [001] distortion (P4mTO) and 150 
meV/f.u. with r^[110] {Amm,2). The PAram phase and 
Amm2 phase are further favored by compressive and ten- 
sile strain, respectively, and their energies drop below 
that of the FM PA/nimm phase at sufficiently large val- 
ues of the strain. In addition to lowering the energy, polar 
distortion of the metallic G-AFM PA/mmm phase tends 
to open a band gap, with the polar (P4mm) phase hav- 
ing a nonzero gap for all strains considered, and Amm2 
having a nonzero gap for tensile strain above -1-1.7%. 

Rotational and JT distortions can also lower the total 
energies of certain phases with G-AFM ordering. How- 
ever, we find that structures with only rotational {Pbnm, 
P2i/m) or only M-type JT distortions (Pi/mbm) of 
the high-symmetry PA/mmm phase are not energetically 
competitive with the cubic FM phase. On the other 
hand, rotational instabilities in the polar AFM PAram 
phase further lower the energy and lower the symmetry to 
Pna2i, lowering the critical strain of the phase boundary 
between FM and AFM from -1.5% to -0.7%. For AFM 
A'mm2, the rotational instabilities lead to a smaller en- 
ergy lowering and a symmetry lowering to Pmc2i, with 
a small reduction in the critical strain. Finally, we com- 
puted the energies of phases with A- or G- type orderings, 
which appear in the computed epitaxial strain phase di- 
agram of SrMnOaQ. In SrCoOa, at each strain value 
these magnetic orderings prove to be higher in energy 
than the FM or G-AFM phases. 

Considering the lowest energy phase at each strain, we 
find three distinct phases in the epitaxial strain phase 
diagram of SrCoOs, with phase boundaries at -0.72% 
and -1-2.04% strain. Next, we consider the computed 
properties of the epitaxially strained phases. In Fig. [2j 
we present the FE polarization, c-lattice parameter, and 
band gap of the lowest energy phase at strains from - 
1.9% to -|-3.3%. The polarizations of the FE phases 
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FIG. 2: (Color online) (a) Computed GGA+U ferroelectric 
polarization of SrCoOs G-AFM (square) in the lowest energy 
structure at each strain value. Open and solid symbols repre- 
sent structures where ferroelectric polarization is along [110] 
and [001] respectively, (b) Band gap (c) c-lattice parameter 
for the lowest energy structure at each strain value. Symbols 
as in Fig. [T] 



are comparable to that of bulk BaTiOs, and increase 
with increasing strain as expected from the polarization- 
strain coupling characteristic of ferroelectric perovskites. 
The large jumps in c-lattice parameter from the EMM 
phase to the FE^ phase can also be attributed to the 
strain-polarization coupling. The computed band gaps 
are above 0.25 eV, and these DET values are expected 
to be substantial underestimates of the yet-to-be mea- 
sured values. The role of the polar distortion in opening 
a band gap was previously noted in the related compound 
SrMnOs [7|. While the band gap is mainly controlled by 
the amplitude of the polar distortion in SrCoOs as well, 
the rotational distortions and strain also affect the value, 
accounting for the epitaxial strain dependence of the gap 
in the FE^ phase. 

Both of the phase boundaries in the computed epitax- 
ial strain phase diagram of SrCoOa are first order, with 
large jumps in properties across the boundary. Thus, for 
a system in the vicinity of one of the phase boundaries, 
large mixed magnetic-electric-elastic responses are pre- 
dicted as the result of the possibility of switching between 
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two phases with quite different magnetic ordering, polar- 
ization, band gap and c-lattice parameter. In particular, 
a AFM-FE phase just at the phase boundary could be 
driven by applied magnetic field to the cubic FM phase, 
with a downward jump in c-lattice parameter; alterna- 
tive, with applied uniaxial compressive stress, a nonzero 
magnetization could be induced. Conversely, a cubic FM 
phase just at the phase boundary could be driven by an 
appropriately oriented applied electric field to the adja- 
cent AFM-FE phase, with a jump in magnetization and 
c-latticc parameter. 

The most interesting behavior associated with these 
phase boundaries is the possibility of a metal-insulator 
transition triggered by applied electric or magnetic fields 
or uniaxial stress. In applied magnetic field, a AFM-FE 
system just at the phase boundary will not only show 
a jump in magnetization, polarization and c-lattice pa- 
rameter, but also a sharp transition from an insulating 
to a metallic state. The jumps in c-lattice parameter 
and polarization also offer the possibility of triggering 
an insulator-metal transition in the AFM-FE state and 
a metal-insulator transition in the metallic FM state by 
applying uniaxial stress and an electric field, respectively. 
The on-off ratio is expected to be large, since the insu- 
lating phase has a band gap of at least several tenths 
of an eV and thus a low conductivity, while the metallic 
phase should have high conductivity since SrCoOa has no 
substitutional disorder that would increase scattering. 

The topology of the phase diagram is not sensitive 
to changes in U. Quantitatively, in the extreme case 
U=0, the computed critical strains are larger (about 
±4%) because the energy difference between FM and 
G-AFM {AE=300 meV) is larger than in GGA+U 
(A£'=:258 mcV), and the frequency shift of Slater phonon 
(w(FM)=163cm-i to u;(G-AFM)=124icm-i) is reduced 
compared to the GGA+U value reported in Table HI 

In experimental investigation of these predictions, 
there will be two major challenges. The first, common to 
all systems with first-order phase transitions, will be to 
minimize the hysteresis in switching from one phase to 
another; one strategy to increase reversibility has been 
demonstrated in Ref. [s^- The other, specific to SrCoOs, 
will be to reduce the oxygen vacancy concentration to 
achieve stoichiometric SrCoOa epitaxial films. However, 
the unusually small elastic energy associated with the 
AFM-FE phases at and above the critical strains should 
greatly facilitate the growth of the strained films. 

In summary, we have presented first-principles calcula- 
tions of the epitaxial-strain phase diagram of perovskite 
SrCoOa. Through combination of the large spin-phonon 
coupling with polarization-strain coupling and coupling 
of the band gap to the polar distortion, both tensile and 
compressive epitaxial strain are seen to drive the bulk 
FM-M phase to AFM-I-FE phases. At these coupled 
magnetic-ferroelectric metal-insulator phase boundaries, 
cross responses to applied electric and magnetic fields and 



stresses are expected. An applied magnetic field can in- 
duce a transition from the AFM phases to the FM phase, 
with accompanying disappearance of polarization, down- 
ward jump in c lattice parameter, and an insulator-metal 
transition, the latter producing a highly nonlinear mag- 
netoresistance. Similarly, an applied uniaxial compres- 
sive stress will change the magnetic ordering and induce 
an insulator-metal transition by driving the AFM-I-FE 
phases to the FM-M phase. Finally, an applied electric 
field in the FM-M cubic phase will drive the system to the 
AFM-I-FE phases, reducing the magnetization to zero 
and inducing a nonzero polarization and metal-insulator 
transition. 
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